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Abstract 

The charmed strange baryon S° was searched for in the decay channel S° — > 
2~7r + , and the beauty strange baryon in the inclusive channel — > 5~£~z/X, 
using the 3.5 million hadronic Z events collected by the DELPHI experiment in 
the years 1992-1995. The S~ was reconstructed through the decay 2~ — > Att~, 
using a constrained fit method for cascade decays. An iterative discriminant 
analysis was used for the S° and selection. The production rates were mea- 
sured to be / H oxBR(H° -> = (4.7 ± lA(stat.) ± l.l(syst.)) x 1(T 4 per 
hadronic Z decay, and BR(6 -> S 6 )xBR(S 6 -> S-M) = (3.0 ± l.O(stoi) ± 
0.3(syst.)) x 10~ 4 for each lepton species (electron or muon). The lifetime 
of the S 6 baryon was measured to be r^ b = 1. 45^43 (stat.) ± 0.13(syst.) ps. 
A combination with the previous DELPHI lifetime measurement gives r S() = 
1.48+°; 4 °(stat.) ± 0.12(syst.) ps. 
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1 Introduction 

Measuring the production rates of baryons, and heavy baryons in particular, is im- 
portant in order to understand the underlying fragmentation process in Z — > qq events. 
The fragmentation process involves small momentum transfers and perturbation theory is 
not applicable, consequently no good theoretical description exists and phenomenological 
models have to be used. The production of a baryon-pair requires the creation of a 
di-quark pair in the fragmentation. The exact nature of the mechanism by which this 
occurs is still largely unknown. Thus the tuning of fragmentation models and the un- 
derstanding of the processes involved in baryon production, require good measurements 
of the production of all the baryons in general, and baryons containing heavy quarks in 
particular. 

In this paper, a first measurement of the production at the Z resonance of the 
charm strange baryon S° is presented 1 , using the exclusive decay channel S° — > S~7r + . 
As a cross-check, the S(1530)° resonance is reconstructed, through the decay channel 
S(1530)° -> 2-7T+. 

A measurement of the production and lifetime of the strange 6-baryon is also 
presented, using the semileptonic decay channel, S& — > S^i^uX. Here is used as a 
notation for the strange 6-baryon states and S°. The baryon will decay to X c Xi~v 
followed by X c — > E~X , where X c is a charmed baryon which yields a H~ hyperon. X c 
is dominantly S°, thus the most common state in this decay channel is the E^. 

A first observation of the St baryon production and lifetime has been published by 
DELPHI, using a smaller data sample [1], and has been confirmed by ALEPH [2]. Here the 
full LEP1 data sample collected by the DELPHI experiment between the years 1992-1995 
is used. In this paper, the background is determined from the data whereas simulation 
was used in the previous DELPHI analysis. The new lifetime measurement is statistically 
independent from the previous DELPHI measurement and relies on a different method. 

A constrained multivertex fit has been performed to reconstruct the S - — > A7r~decay. 
For the S°, 5(1530)° and E b selections, an iterative discriminant analysis method has 
been applied. 

2 The apparatus 

The DELPHI detector is described in detail elsewhere [3,4]. The detectors most im- 
portant for this analysis are the Vertex Detector (VD), the Inner Detector (ID), the Time 
Projection Chamber (TPC), and the Outer Detector (OD). For the lepton identification 
in the analysis the electromagnetic calorimeter (HPC) and the muon chambers were 
also used. 

The VD consists of three concentric layers of silicon-strip detectors, located at radii 
of 6 cm, 9 cm and 11 cm from the beam axis. The polar angles 2 covered for particles 
crossing all three layers are 44° < 9 < 136°. In 1994 and 1995, the first and third layers of 
the VD had double-sided readout and gave both Exf> and z coordinates. The TPC is the 
main tracking device where charged-particle tracks are reconstructed in three dimensions 
for radii between 40 cm and 110 cm. The ID and OD are two drift chambers located 
at radii between 12 cm and 28 cm and between 198 cm and 206 cm, respectively, and 
provide additional points for the track reconstruction. 

1 Charge conjugated states arc implied throughout this paper, unless otherwise stated. 

2 In the standard DELPHI coordinate system, the z axis is along the electron beam direction, the x axis points towards 
the center of LEP, and the y axis points upwards. T he polar angle to the z axis is called and the azimuthal angle around 
the z axis is called <j>; the radial coordinate is R = \J x 2 + y 2 . 
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The 6-tagging package developed by the DELPHI collaboration [4,5] has been used to 
select Z — > bb events. The impact parameters of the charged-particle tracks, with respect 
to the primary vertex, have been used to build the probability that all tracks come from 
the primary vertex. Due to the long 6-hadron lifetime, the probability distribution is 
peaked at zero for events containing 6-quarks whereas it is flat for events containing light 
quarks only. 

3 Event selection and simulation 

Hadronic Z decays were selected by requiring at least four reconstructed charged par- 
ticles and a total energy of these particles (assumed to be pions) larger than 12% of the 
centre-of-mass energy. The charged-particle tracks had to be longer than 30 cm in the 
TPC, with a momentum larger than 400 MeV/c and a polar angle between 20° and 160°. 
The polar angle of the thrust axis, Q thrust, was computed for each event and events were 
rejected if | cos 8thrust\ was greater than 0.95. With these requirements the efficiency for 
the hadronic Z selection was larger than 95%. A total of 3.5 million hadronic events were 
selected. 

Simulated events were produced by the Jetset 7.3 parton-shower generator [6] and 
then processed through a detailed simulation program, Delsim, which modelled the 
detector response [4] . The simulated result from Delsim was then processed by the same 
reconstruction program as used for the data, DELANA [4]. A total of 9.8 million Z^ qq 
events was simulated (11.8 million for the E b analysis). 

In this simulation sample, there were only a few thousand S° — > S~7r + events, and 
about 1000 E b — > EiuX events. Thus some dedicated S° and E b samples were generated 
for the years 1992-1995 using the Delsim and Delana versions corresponding to each 
year. About 58 000 H° — > E^tt + events and 208 000 S& — > EiuX events were simulated 
(see Table 1). 



Year 


S" events 


E b events 


Sjj events 


1992 


11 159 


50 710 


1 365 


1993 


12 142 


52 735 


1 331 


1994 


20 636 


50 203 


1 221 


1995 


13 874 


49 378 


1 227 


Total 


57 811 


203 026 


5 144 



Table 1: Number of simulated S° — * S ix + , E b and E b events for each 1992-1995 year. 



4 S reconstruction 

The S~ hyperon was reconstructed through the decay S _ — > Att~. A constrained 
multivertex fit to the three-dimensional decay topology was used to reconstruct the decay 
chain and suppress the combinatorial background. 

In this analysis all V° candidates, i.e. all pairs of oppositely charged particles, were 
considered as A candidates. For each pair, the highest momentum particle was assumed to 
be a proton and the other a pion, and a vertex fit was performed by the standard DELPHI 
V° search algorithm [4]. The A candidates were selected by requiring: an invariant mass 
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M{pn~) between 1.107 GeV/c 2 and 1.125 GeV/c 2 , a \ 2 probability of the V° vertex 
fit larger than 0.001 and an R(p decay length greater than 1.0 cm. To avoid gamma 
conversions, the relative transverse momentum, pt, of the proton and pion had to be 
greater than 0.03 GeV/c, with pr calculated with respect to the line joining the primary 
and secondary vertices. The angle in the rcy-plane between the V° momentum and its 
line of flight had to be smaller than 0.08 radian. If the V° was reconstructed outside 
the VD, it was also required that no signal in the vertex detector could be consistently 
associated with the V° vertex tracks. In the following, the proton and pion from the A 
decay will be called p\ and 7Ti, respectively, thus denoting A = (p^^). 

The A's selected as described above were combined with pions (called 7r 2 in the follow- 
ing) with the same charge as the 7Ti from the A. The 7r 2 candidate had not to be tagged 
as an electron or a muon. 

A constrained multivertex fit [7,8] was performed if the invariant mass M{Kt^2 ) was 
smaller than 2.0 GeV/c 2 and if the distance between the two trajectories in the z direction 
was smaller than 2 cm at the point of crossing in the xy-plane. The fit used was a general 
least-squares fit with the following kinematical and geometrical constraints applied to 
each S~ candidate: 

- the invariant mass M(pi7rj~) had to be equal to the nominal mass of the A; 

- the Rep and z coordinates of p\ and i\\ had to be the same at the radial distance of 
the A decay point; 

- the R(f> and z coordinates of A and tt 2 had to be equal at the radial distance of the 
2~ decay point; 

- to ensure momentum conservation at the E~ decay point, the polar and azimuthal 
angles of the S~ candidate had to be equal to the angles from the curved trajectory 
between the decay and primary vertices. The curvature was calculated from the S~ 
momentum. 

For each of the three particles, i.e. the proton and pion from the A and the pion from 
the S~, the fit was performed with the following five track parameters: 1/r (r being the 
radius of curvature of the track), the z and R(f> impact parameters, the polar angle 9, and 
the azimuthal angle 0. These 15 variables, plus the z coordinate of the primary vertex, 
were the measured variables in the fit. The x and y coordinates of the primary vertex 
were so precisely measured that they were taken as fixed. The unmeasured variables were 
the radial distances of the S" and A decay points, giving a total of 18 variables in the 
fit. The curved S~ track was not measured, but calculated in the fit. All the tracks were 
corrected for ionization losses, according to the given mass hypothesis. The performance 
of the fit was tuned by adjusting the covariance matrices of the tracks. The adjustment 
consisted in a scaling of the errors of the track parameters. After the adjustment the pull 
distributions of the 16 fitted quantities were standard normal distributions within 10%. 

The events for which the fit converged gave the Air^ invariant-mass spectrum shown in 
Figure 1, The solid line is a fit, for illustrative purpose, using two Gaussian functions with 
the same mean for the signal, and a first order polynomial for the background, yielding 
9445±584 S~. 

5 Iterative discriminant analysis 

Using sequential cuts to select the signal leads to what can be pictured as a multi- 
dimensional rectangular box in the parameter space. A more flexible separating surface 
can be obtained if the variables are combined in a polynomial instead. In this analysis a 



4 



non-linear discriminant D of the form 

D = x T (a + Bx) = aixi + ... + a n x n + B lx x\ + ... + B Xn x x x n + ... + B nn x 2 n (1) 

was defined [9], where Xi are the n variables used, and ai and are the corresponding 
weights. This can be written as 

D = c T y = ciyi + ... + c m y m , (2) 

where c contains all the weights and -B^, y is the vector (x, x ■ x T ), and m = (n 2 + 
3n)/2. To obtain a good signal- from-background separation, the variance of D should be 
as small as possible while the separation in D should be as large as possible. Therefore 
the ratio 

was maximized. Here V is the sum of the signal and background variance matrices 

V = v.,, + V bkg (4) 

and A/x is the difference between the signal and background arithmetic means: 

A/x = V sig - » bkg (5) 

of all the m variables in y. The discriminant was calculated iteratively, and at each 
step of the iteration the variables Xi that gave the maximum background rejection were 
chosen. A chosen signal efficiency was used to determine the D-cuts after each step. The 
total number of variables used in the discriminant D, as well as the number of iterative 
steps, and the efficiency of each step, can be varied in the discriminant analysis. A 
large number of variations of the combinations of variables used in the discriminant, the 
number of steps, as well as the signal efficiencies, was investigated before deciding on 
which parameter combination to use. 



6 Production of S tt + states 

Each S~ candidate in the mass range 1.30 < M(An 2 ) < 1.34 GeV/c 2 , was combined 
with another pion, called 7r 3 . It was required that 7r 3 should have a charge opposite to 
the tti from A decay, and should not be tagged as a lepton. For all (S~7T3 ) combinations 
the invariant mass was calculated. If its value satisfied 2.2 < M(S~7r 3 f ) < 2.75 GeV/c 2 
(for S° ), or M(S _ 7T3 ) < 1.6 GeV/c 2 (for S(1530)°), the iterative discriminant analysis 
described above was performed. 

6.1 S[! — > S 7r + selection 

6.1.1 Separate simulation sets 

Two separate sets of simulated events were needed. The training sample (called T(s°) 
in the following) was used to find the weight vector c that gave a maximum signal- 
from-background separation in D, while the analysis sample (called A(h 1 c »)) was used to 
determine the S° efficiency (see Section 6.1.3). Both sets T(h°) and A(s c c ') consisted of signal 
as well as background events. The dedicated S° — > H~7r + events described in Section 3 
were thus divided into two sets: about 2/3 of the events were used in set T(h°) and the 
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rest in set A(s°). The simulated Z— > qq events were used for the background in both 
cases, using only events without any simulated S°'s. The signal in set T(h°) consisted 
of about 36000 events, while the background corresponded to approximately 4.4 ■ 10 6 
qq events. After the constrained fit, and the above-mentioned selections, about 1900 
signal and 75000 background events remained in set T(s°) to be used for the discriminant 
training. Set A(H°) consisted of approximately 22000 2° — > events, plus background 

events without any simulated S° corresponding to 5.3 ■ 10 6 qq events. For each year, the 
number of background events in set A(3°) was weighted to correspond to the number of 
data events. Each of the years 1992-1995 was trained separately. The total numbers of 
events in the samples used for the (S°) analysis are shown in Table 2. 





MC T(~<>) 


MC A(h<>) 




Data 


Signal 
After cuts 


35 678 
1 883 


22 133 
1 268 






Background 
After cuts 


4 426 483 
74 750 


5 329 144 
86 967 


Total 
After cuts 


3 498 492 
57 291 



Table 2: The numbers of signal and background events in the simulation sets T"(h°) and 
A(s°) (described in Section 6.1.1), and in the data, both before and after the multivertex 
fit and applied selections. 



6.1.2 Parameter optimization 

As mentioned above, a large number of parameter combinations was investigated be- 
fore a discriminant with two iterative steps, and seven variables, was chosen. The chosen 
variables were the momentum p(S°), the fitted mass M(E~)f it , the angle between S~ and 
7i\j~, the 6-tagging probability, the A decay length or the S~ decay length, the momentum 
p(S~) or the A mass, the impact parameter of the nT, or the momentum p(7rj), with 
some variations among the different years. 

6.1.3 E° c production 

Using the discriminants D with the above mentioned variables on the simulation set 
A(s°), resulted in the (An^n + ) invariant-mass distribution shown in Figure 2. An ex- 
tended unbinned maximum-likelihood fit was performed. The likelihood function used 
had the form: 

log L = log (jJL a f a + fibfb) - {(i s + fi b ) , (6) 

where f s is a Gaussian function used to parametrize the signal and /*,, the probability 
density function for the background, is a first order polynomial. The numbers of signal 
and background events were given by fi s and fib, respectively. The fit gave 498 ± 28 
events with a a width of 19±2 MeV/c 2 and mass M(S°) = 2471 ± 1 MeV/c 2 , which can 
be compared with the generated mass of 2473.0 MeV/c 2 . The number of true S° events 
in the peak was 494, and there was thus an excellent agreement. 

The same discriminants were also used on the data, with the resulting mass distribution 
shown in Figure 3. In this case the fit gave 45±13 events with a a width of 22±6 MeV/c 2 . 
The mass from the fit was 2460±8 MeV/c 2 , in agreement with the world average value 
M(S°) = 2471.8 ± 1.4 MeV/c 2 [10]. 
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The reconstruction efficiency is highly momentum dependent, and in order to avoid 
biases due to the S° momentum distribution not having been correctly described by the 
simulation, the S° rate was determined in two momentum bins, as shown in Table 3 and 
Figure 4. Integrating the observed distribution and using the JETSET generator [6] to 
estimate the fraction of events outside the observable momentum region (5.5%), a total 
rate of 

/ H o x BR(S° -> S-7T+) = (4.7 ± lA(stat.)) x 1(T 4 (7) 
per hadronic Z decay was obtained. 



Particle 


x p interval 


Efficiency 


Nb. of particles 


1 dN 

N h ad dx p 




0.10-0.40 
0.40-0.80 


(2.9±0.2)% 
(1.9±0.2)% 


28.2±10.6 
11.0±6.5 


(9.2 ±3.5) x 10" 4 
(4.1 ±2.4) x 10" 4 


H(1530) u -> 2-7T+ 


0.015-0.10 
0.10-0.20 
0.20-0.50 


(3.4±0.5)% 
(7.6±0.9)% 
(1.5±0.4)% 


271.1±39.3 
293.1±37.5 
37.3±15.8 


(2.7 ±0.4) x 10~ 2 
(1.1 ±0.1) x 10~ 2 
(2.4 ± 1.0) x 10~ 3 



Table 3: The H° and S(1530)° differential production rates. The efficiency and number 
of particles from the fit in each of the x p = p/pbeam intervals are also given. 



6.2 S(1530)° production 

As for the S° analysis, an iterative discriminant analysis was applied for the S(1530)° 
selection, and two separate simulation sets were used. No special S(1530)° simulation was 
needed, and the qq simulated events were used for both signal and background events. 

The S(1530)° discriminant analysis was applied to the simulation and resulted in the 
(A7r~7r + ) invariant-mass distribution shown in Figure 5. The same maximum likelihood 
function as for the S° was used, except that the signal was described by a Breit-Wigner 
function, while the background was parametrized by [11] 

F(x) = (x - x ) a x exp (b (x - x ) +h(x- x ) 2 ), (8) 

where x is the invariant mass measured in GeV/c 2 , xo is the kinematical limit of 
1.4609 GeV/c 2 , and a, bo and b\ are free parameters. Since the signal peak is fairly close 
to the maximum of the background, the fit easily became unstable. Therefore the width 
T in the Breit-Wigner was kept fixed at the world average value of 9.1±0.5 MeV/c 2 [10]. 

Performing the fit on the simulation resulted in 844 ± 71 signal events, and a mass 
of 1531.2 ± 0.5 MeV/c 2 , in good agreement with the number of S(1530)° events in the 
sample (831), and the generated mass of 1532.0 MeV/c 2 . 

Using the same discriminant on the 1992-1995 data sample gave the invariant mass 
spectrum shown in Figure 6. The unbinned maximum-likelihood fit gave 599 ± 57 signal 
events in the peak, with a mass of 1533.0 ± 0.5 MeV/c 2 , which can be compared with the 
world average value M(S(1530)°) = 1531.80±0.32 MeV/c 2 [10]. The S( 1530)° production 
rate was evaluated in three momentum bins, as shown in Table 3 and Figure 4. Integrating 
the observed distribution and using the JETSET generator [6] to estimate the fraction 
of events outside the observable momentum region (7.7%), a total rate of 

/s(i530)° x BR(S(1530)° -> S~7r + ) = (4.5 ± 0.5) x 10" 3 (9) 

per hadronic Z decay was obtained. 



7 



6.3 Systematic uncertainties in production fractions 
6.3.1 The 2° baryon 

Several different sources of systematic uncertainties were investigated. The systematic 
uncertainty from the choice of the discriminant parameters was studied by varying the 
number of discriminant variables, the number of steps, and the signal efficiencies. No 
significant variation beyond the expected statistical fluctuation was found. 

Due to the difference in the S° momentum distribution in c — > S° and b — > c — > S° 
events, respectively, the reconstruction efficiencies in cc and bb events can also differ. 
An important contribution to the systematic error was therefore the uncertainty in the 
different c and b efficiencies, as well as the relative production of S° in cc and bb events. 
The number of observed S° particles, N seen , is given by 

N seen = N X (R c • f c • € c + R b ■ f b ■ 6 b ) X BR, (10) 

where N is the total number of data events, R c and R b are the Z partial widths 
r(cc) /T(hadrons) and T(bb) /T(hadrons) , e c and e b are the S° reconstruction efficien- 
cies for cc and bb events, f c = f(cc — > S°) and f b — f(bb — > S°) are the fractions of cc and 
66 quark pairs giving a S°, respectively (it has been assumed that S° is only produced 
in cc and bb events). Finally, BR = BR(S° — > H _ 7r + ) is the branching fraction for the 
studied channel. 

The reconstruction efficiencies for cc and bb events were determined from the simulation 
and found to be e c = (1.4±0.2)%, and e b = (3.2±0.2)%. To estimate the relative weights 
of the f c and f b terms, both simulation and data events were used. The 6-purity of the 
S° sample was enhanced with a 6-tag probability selection cut and f b x BR could then be 
found assuming that the contribution from cc events was negligible. The new b efficiency 
was found from the simulation. The number of fitted S° in data after this selection, 
Nseln, was measured, and used to determine f b x BR from: 

Nill = Nx(R b -f b -e ( b b) )xBR. (11) 

It was found that f b x BR = (1.3 ± 0.6) x 10~ 3 , essentially independent of the chosen 
6-tag cut. 

Using this value in equation 10, j c j f b was found to be compatible with 1. For the study 
of the systematic uncertainty f b x BR was varied by one statistical standard deviation, 
while f c / fb was varied between 0.5 and 2.0 in equation 10. The systematic uncertainty 
was thus estimated to be ±1.0 x 10~ 4 . 

Other sources of systematics came from the contribution of the finite simulation statis- 
tics to the uncertainty on the total S° efficiency, and the error rescaling done in the 
multivertex fit, giving a systematic uncertainty of ±0.3 x 10 -4 from each source. 

As already mentioned, the JETSET model was used to estimate the fraction of events 
in the unobserved momentum regions, which was found to be 5.5%. The comparison 
with JETSET (Figure 4) shows a good agreement, thus the systematic uncertainty was 
estimated by varying the value from the simulation by ±50%, resulting in a systematic 
uncertainty contribution of ±0.2 x 10~ 4 . 

All the above systematics are summarized in Table 4. These uncertainties were then 
added in quadrature, which gave the final result: 



/ s o x BR(S° -> = (4.7 ± lA(stat.) ± l.l(syst)) x 10~ 4 . 



(12) 
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Source 


/xBR(S^) 


/xBR(S(1530)°) 


Finite MC statistics 
MC extrapolation 
multivertex fit 
b, c efficiencies 


0.3xl0~ 4 
0.2xl0" 4 
0.3xl0" 4 
l.OxlO" 4 


0.4xl0" 3 
0.2xl0- 3 
0.4xl0- 3 


Total 


l.lxlO^ 4 


0.6xl0~ 3 



Table 4: The different contributions to the total systematic uncertainty for S° and 
S(1530)° 7 as described in Section 6.3. 

6.3.2 The 5(1530)° baryon 

The systematic uncertainties for H(1530)° were studied in the same way as for the 
S° events, except that since the S(1530)° mainly comes from the fragmentation, these 
results were not affected by the flavour of the leading quark. All the different uncertainty 
contributions, summarized in Table 4, were then added in quadrature, and the final 
2(1530)° result was 

/s(i530)° x BR(S(1530)° -> 2-7T+) = (4.5 ± 0.5(stat.) ± 0.6(syst.)) x 10" 3 , (13) 

in agreement with the world average value of (5.3 ±1.3) x 10~ 3 [10] and the results 
of DELPHI [12] and OPAL [11]. The 2(1530)° result was used as a cross-check of the 
analysis method. 

7 Update of the E production rate and lifetime 

The strange 6-baryon H5 was searched for in the semileptonic decay channel, — > 
S^i^uX. In the semileptonic decays of heavy hadrons the flavour of the spectator system 
of the initial state is transmitted to the final state. This property can be used to study 
Sf, baryons from the observation of S T production accompanied by a lepton of the same 
sign. The occurrence of S T — £ T pairs of same sign ("right sign") is then compared to 
that of opposite sign pairs, S T — £ ± ("wrong sign"). 

7.1 H~ and lepton reconstruction 

The S~ was reconstructed using a constrained multivertex fit as in the S° analysis. If 
the fit was successful, the H~ candidate was combined with a lepton candidate (electron 
or muon) within 1.0 radian of the H~ momentum vector. Since the expected production 
rate of the E b is very small, loose selections were applied to the H~ and A candidates. 
The discriminant analysis method described above was used for the final E b selection. 
Five variables were used in the discriminant; the transverse momentum of the lepton 
with respect to the jet axis, the invariant mass of the S~ and lepton, the combined 
momentum of the S~ and lepton, the number of charged particles in a 0.31 radian cone 
around the lepton direction and the S~ variable, £ = — lnx p , where x p = p-=/Vbeam,- 

When applied to the Monte Carlo analysis sample, consisting of 1/3 of the simulated 
events, the discriminant method gave the resulting An invariant-mass distributions of Fig- 
ure 7. Applying the extended unbinned maximum-likelihood fit described in Section 6.1.3 
to the two distributions, with a Gaussian function for the S~-peak and a constant value 
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for the background, gave 34.2±5.9 right-sign events, and 11.3±3.5 wrong-sign events, 
when normalised to the size of the data sample. The number of true E b events in the 
right-sign E~ mass peak was 25.6. Here the mass peak region was defined as the region 
with reconstructed E~ mass within ±10 MeV/c 2 of the nominal E~ mass. 

The same analysis was applied to the full DELPHI 1992-1995 data sample and 
the resulting Air invariant-mass distributions are shown in Figure 8. The unbinned 
maximum-likelihood fit to the two distributions, gave a mean value of 1321. 0±0. 8 MeV/c 2 
for the right-sign distribution, compatible with the nominal E~ mass (1321.31 ± 
0.13 MeV/c 2 [10]). The mean value from the fit to the right-sign distribution was used 
as a fixed parameter in the fit to the wrong-sign distribution. The maximum-likelihood 
fit resulted in 28.3±5.8 right-sign events, and 7.6±3.3 wrong-sign events. 

7.2 Eb production rate 

The number of background events in the right-sign E~ mass peak was estimated from 
that in the wrong-sign mass peak, which resulted in 20.7 ± 6.7 E b events found in the 
data. According to the simulation, the number of background events is equal in the right- 
and wrong-sign E~ mass peaks. In the Monte Carlo analysis sample the same procedure 
gave 23.0 ± 6.8 events to be compared with the number of true E b — > E~£~X events in 
the right-sign mass peak which was 25.6. 

The total E b efficiency, calculated from the E b signal simulation sample, was 
(2.3 ± 0.1)%. Using the measured fraction of bb relative to all Z hadronic decays, R b 
= (21.650 ± 0.072)% [10], leads to a E b production rate of: 

BR(b -> E b ) x BR(E b -> 5TX) = (3.0 ± 1.0(stat.)) x 10~ 4 (14) 

per lepton species, averaged for electrons and muons. 

The dominating source of systematic uncertainty was the contribution of A b to the 
background. The uncertainty in the background contribution of A b was estimated by 
varying the amount of these events in the background by ±20% of the value in the sim- 
ulation [10]. This gave a shift of ±0.20 x 10~ 4 of the production rate. Other sources 
of systematic uncertainty were the finite simulation statistics, the error rescaling done 
in the multivertex fit, the Monte Carlo extrapolation of the events into the unobserved 
momentum regions (9.3%) and a possible E b polarisation. All the above sources of sys- 
tematic uncertainty are summarized in Table 5. The final result for the E b production 
rate was: 

BR(6 -> E b ) x BR(E b -> E~TX) = (3.0 ± l.O(stat) ± 0.3(syst.)) x 10~ 4 (15) 

per lepton species, averaged for electrons and muons. This measurement of the production 
rate is in agreement with the previous measurements done by DELPHI [1] using a smaller 
data sample, and by ALEPH [2], see Table 6. In the previous DELPHI analysis [1], the 
background was estimated from the simulation while in this analysis the background was 
estimated using the wrong-sign data sample. 

7.3 Eb lifetime 

The final sample of E~l~ events is also used to measure the E b lifetime. The secondary 
vertex from the E b semileptonic decay is obtained by use of the BSAURUS package [13]. 
The secondary vertex in the hemisphere is calculated in BSAURUS using tracks that are 
likely to have originated from the decay chain of a weakly decaying 6-hadron state. The 
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Source 


Production rate variation 


Finite MC stat. 


O.lOxlO- 4 


Multivertex fit 


0.15xl0" 4 


MC extrapolation 


0.15xl0" 4 


E b polarisation 


0.14xl0" 4 


MC model of background 


0.20X10" 4 


Total 


0.3xl0" 4 



Table 5: The different contributions to the total systematic uncertainty of the S& produc- 
tion rate, as described in Section 7.2. 



Reference 


Production rate 


ALEPH [2] 
DELPHI [1] 


(5.4 ± l.l(stat) ± O.B(syst.)) x 10~ 4 
(5.9 ± 2.1(stat.) ± 1.0(syst.)) x 10~ 4 


this analysis 


(3.0 ± l.O(stat) ± 0.3(syst.)) x 10" 4 



Table 6: Comparison between different results for the E b production rate. 



vertex fitting is done in three dimensions, using the constraint of the direction of the 
6-hadron candidate momentum vector, p b . The decay-length estimate in the R(j) plane, 
is defined as the distance between the primary- and secondary-vertex positions if the 
secondary-vertex fit was successful. The three-dimensional decay length is obtained as 
L = Lr^/ sin6, where 8 is the polar angle of the 6-hadron candidate. The sign of the 
decay length is determined by the direction of the S£ momentum vector: the distance is 
positive if the secondary vertex is found beyond the primary vertex in this direction. The 
resulting proper-time estimate for the E b candidates is given by t = L ms b /p b , where the 
E b rest mass is taken to be 5.8 GeV/c 2 . 

The lifetime was determined by an unbinned maximum-likelihood fit to the proper 
time distribution. 

The background was taken as the wrong-sign combinations plus the right-sign combi- 
nations outside the mass peak but with a reconstructed S~ mass between 1.280 GeV/c 2 
and 1.363 GeV/c 2 . Both signal and background lifetimes were fitted simultaneously. The 
purity in the right-sign mass peak was fixed in the fit and was taken from the simulation 
to be 0.67. 

In the 1992-1995 DELPHI data sample, 29 events in the right-sign mass peak remained 
after the secondary vertex fit. The statistical overlap of one event with the previous 
DELPHI lifetime measurement was removed. The result of the lifetime fit on the DELPHI 
data is presented in Figure 9, and gave r= 6 = 1.45loH(s£ai.) ps. The exact composition 
of the background, and the lifetimes of its individual components, had no effect on the 
signal lifetime since the background lifetime was fitted on the data. Varying the fraction 
of Sft in the right-sign peak between 0.60 and 0.74 resulted in a variation of the fitted S b 
lifetime of ±0.10 ps. The proper time resolution was varied by ±50%, which resulted in 
a shift of ±0.07 ps. The effect of a possible E b polarisation was studied and found to be 
small. Finally the E b lifetime was measured to be 

r Sb = lA5t° 5 i(stat.) ± 0.13(syst.) ps. (16) 
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The measurement is in agreement with the previous measurements done by DELPHI [1] 
and by ALEPH [2], see Table 7. 



Reference 


Lifetime (ps) 


ALEPH [2] 
DELPHI [1] 


1.35^(^0^(^0 
1.5t° l(stat.)±0.3(syst.) 


this analysis 


1.45+^(stat.) ±0.13(syst.) 



Table 7: Comparison between different results for the S& lifetime. 

The earlier DELPHI lifetime measurement [1] used the data from 1991-1993 and a 
different method to reconstruct the H~-hyperon and the proper time, and the background 
lifetime was estimated using simulation. A combination of the two DELPHI lifetime 
measurements gives 

r Sb = lASt° ii(stat.) ± 0A2(syst.) ps, (17) 
using the method outlined in [14]. The systematics are uncorrected. 

8 Summary and conclusions 

The production rate per hadronic Z decay for the charmed baryon S° has been mea- 
sured for the first time: 

/sg x BR(H° -> S~7r + ) = (4.7 ± lA(stat.) ± l.l(syst.)) x KT 4 . 

As a cross-check, the S(1530)° resonance was also reconstructed, and the corresponding 
production rate was found to be: 

/e(i530)° x BR(S(1530)° -> E~tt + ) = (4.5 ± 0.5(stat.) ± 0.6(syst.)) x 10~ 3 , 

in agreement with previous results [10]. 

The beauty strange baryon E b was searched for in the semileptonic decay channel 
E b — > E~£~X. The product of the branching ratios in b and E b decays was measured to 
be: 

BR(b -> E b ) x BR(E b -> ETX) = (3.0 ± 1.0(stat.) ± 0.3(syst.)) x 10" 4 

per lepton species, averaged for electrons and muons. 
A measurement of the E b lifetime gave: 

r Sb = lA5t° 5 i(stat.)±0A3(syst.) ps, 

in agreement with earlier results [1,2]. A combination of the two DELPHI lifetime mea- 
surements gives 

r Sb = lA8t° Q H(stat.) ± 0A2(syst.) ps. 
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Figure 1: The An~ invariant-mass spectrum, using a constrained fit on the 1992-1995 
data sample. The curve is the result of a fit using a first-order polynomial to parametrize 
the background and two Gaussian functions of same mean for the signal. 
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Figure 2: T/ie An^n + invariant-mass spectrum in the 1992-1995 simulation. The back- 
ground (cross-hatched histogram) corresponds to 5.3 • 10 6 Z— > gg simulated events. The 
signal (white histogram) consists of about 22000 S° — > 2~7r + simulated events. The fit- 
ted curve, as described in Section 6.1.3, uses a first-order polynomial to parametrize the 
background and a Gaussian function for the signal. 
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Figure 3: Same as Figure 2 for the 1992-1995 data sample, corresponding to 3.5 • 10 6 
hadronic Z decays. 
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Figure 4: T/ie measured differential production rates for S° (plotted as triangles), and 
S(1530)° (circles). For comparison the JETSET 7.4 [6] prediction is shown, with a full 
line for and a dashed line for S(1530)°. 
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Figure 5: The An^n + invariant-mass spectrum, using 5.3 • 10 6 simulated hadronic Z 
decays. The background events are cross-hatched. The fitted curve, as described in Sec- 
tion 6.2, uses a parametrized background and a Breit-Wigner function for the signal. 
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Figure 6: Same as Figure 5 for the 1992-1995 data sample, corresponding to 3.5 • 10 6 
hadronic Z decays. 
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Figure 7: T/ie A7r~ invariant-mass spectrum, using 3.8 million simulated hadronic Z 
decays: a) E T — £ T right-sign pairs, b) S T — ^ ± wrong-sign pairs. The true simulated 
events are grey hatched. The two fitted curves, as described in Section 7.1, each uses a 
constant value to parametrize the background and a Gaussian function for the E~ peak. 
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Figure 8: Same as Figure 7 for the 1992-1995 data sample, corresponding to 3.5 ■ 10 6 
hadronic Z decays. 
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Figure 9: The result of the lifetime fit to the selected S^, events in the data sample. The 
dotted curve is for the background, and the dashed line corresponds to the signal. The full 
line is the total. 



